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We study ways of implementing Fermi momentum in elastic 
vector-meson production, and find that the usual on-shell assump- 
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1 Introduction 



The data from HERA have now reached a remarkable level of accuracy, and in 
particular they provide us with excellent information on vector meson elastic 
production [0], 0- We have precise measurements of the photon Q 2 , t and w 2 
dependences of the cross section, for different mesons, as well as information on 
its transverse and longitudinal components, Q 2 being the absolute value of the 
photon off-shellness, t the square of the momentum transferred to the proton and 
w 2 the square of the center-of-mass energy of the 7 ( -*- ) p process. 

Theoretically, this is related to a measurement of the off-diagonal deep-in- 
elastic tensor W^ v {x, Q 2 , m v , t), for t non zero. One may hope that the presence 
of two new scales - the vector meson mass my and the momentum transfer t 
would enable one to concentrate on regions where the perturbative calculation 
is reliable, and observe how one can extend it to smaller scales. Theoretically, 
these studies have recently been set on a firm ground as Collins, Frankfurt and 
Strikman || have shown that a factorisation theorem exists for exclusive vector 
meson production. This theorem asserts that, to leading power of Q 2 and to all 
logarithms, the cross section is the convolution in longitudinal momentum of a 
hard scattering amplitude, an off-diagonal structure function and a meson wave 
function. Furthermore, the theorem can be proven irrespectively of the meson 
mass. Hence at high Q 2 , even light meson production can be reliably calculated. 

The w 2 dependence of the cross section then results from the x, t/w 2 and 
rriy/w 2 dependence of these new distribution functions. It is to be noted that 
these are new objects which are different from the diagonal structure functions. 
Also, the theorem naturally holds for longitudinal vector meson production, which 
is leading in Q 2 , and for which one can show that the off-shell distributions are 
related to the usual gluon distribution. The transverse cross section o~t is a 'priori 
down by a factor Q 2 . However, the factorisation theorem envisions the possibility 
of a different behaviour, but the price to pay, in the case of massless quarks, is 
to allow the leading behaviour of <jt to come totally from the non-perturbative 
region, for quark off-shellnesses of order my /Q 2 . In this paper, we shall go beyond 
this analysis by including massive quarks in a model describing the transition from 
a photon to a QCD bound state. One indeed needs to disentangle the effect of this 
transition on the transverse cross section, before one can extract the off-diagonal 
distributions accurately. 

One expects (and we confirm) that the transition 7* — > V will not modify the 
w 2 dependence significantly at high energy. The w 2 dependence is surely one of 
the outstanding problems of QCD, and our theoretical understanding of it seems 
to have made a big step recently || , although it is difficult to assert yet in which 
direction. We shall make the simplifying assumption that it enters as a constant 
factor which does not depend on any variable but t. We shall come back later to 
the significance of this assumption. 

We could of course, as often done nowadays, assume that this factor is related 
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to W(t = 0) and is proportional to [xg(a;)] 2 . However, we want first to point out 
that this correspondence exists only for longitudinal cross sections, and is only 
approximate. Indeed, as we shall explain below, the imbalance between initial and 
final states in vector-meson elastic production brings one to a different kinematic 
domain from DIS. This has important consequences for quantities related to the 
upper loop of a ladder, especially in the transverse case. 

Hence we shall simply calculate up to a factor, and we shall not claim that we 
can reliably predict it. On the other hand, the Q 2 and my dependences mainly 
(up to logarithms) come from the upper loop. As we shall see, the data is not 
precise enough (or Q 2 not large enough) to require the inclusion of evolution 
effects. As for the t dependence, it comes from the proton form factor, the 
7*V loop, and W. We shall show here that simple assumptions enable one to 
reproduce the data at small t, and that high-t points do not seem to require either 
a pomeron slope or a BFKL enhancement. 

Our starting point is a very simple model that we published recently ||. It 
is based on lowest-order perturbative QCD and on a very naive approximation 
to the meson wave function ||, where in its rest frame the meson is simply 
described by two quarks at rest. Such an approximation should be valid provided 
that the meson motion is negligible, i.e. at very high Q 2 . Despite its simplicity, 
this model was surprisingly successful and reproduced most of the features of 
the data: the mass dependence of the cross section came out naturally, as well 
as the Q 2 dependence. The model also predicted correctly s-channel helicity 
conservation, and even worked in photoproduction. It provided a test of two 
ideas: first of all, elastic production is mainly due to two gluons interacting with 
the meson, secondly, as the mesons were supposed to be made of their constituent 
quarks, the main production mechanism involves only the lowest Fock state. 

Despite its many successes, this model failed, as most other, in the fact that 
the ratio o^/or was predicted to be a factor 8 higher (at Q 2 = 20 GeV 2 ) than 
observed experimentally. As is much larger than cx^, the high-Q 2 total cross 
sections cjl + o T were unaffected by this failure, but it is a strong indication that 
the model is far from complete. We want to emphasize that so far no model 
manages to reproduce the plateau observed experimentally, and that the best 
one has been able to do is to get a slower linear rise. 

We shall show here that a re-consideration of the role of Fermi momentum 
enables one to preserve all the previous successful results while reproducing the 
ratio o"i/o"T. We keep the two fundamental assumptions of our previous paper, 
that the process is dominated by the coupling of two gluons to the lowest Fock 
state. In Section 2, we introduce the general formalism which we shall use for 
the description of the quarks contained in the meson. In sections 3 and 4, we 
give the general form of the amplitude. In section 5, we show that the amplitude 
is in general infrared finite, and in section 6 we give our results, and discuss 
the off-shell contribution to the amplitude. We then conclude, and pose several 
questions concerning the validity of the wave function formalism, and of the 
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connection between diagonal and off-diagonal DIS tensors. 

2 Vertices vs. wave functions 

We shall assume here a form of factorisation, in that we shall not consider soft 
exchanges between the proton and the meson, but only concentrate on a hard 
scattering, a meson vertex function, a singlet exchange and a proton form factor. 
We do keep the exact kinematics, including transverse motion and quark masses, 
in the amplitude describing the transition 7* — > V as we want to assess the 
contribution of near-shell partons. 

We need a model for the off-diagonal structure function, and we adopt the 
simplest one: at high energy, the process is dominated by pomeron exchange, 
which we model by two gluons times a Regge factor. The model of the off diagonal 
structure function will be obtained by convoluting the exchange with a proton 
form factor, which partially kills the infrared divergences. 

The hard process generating the meson, and the qq —>■ V amplitude are treated 
together through the introduction of a meson vertex function. We shall make 
the assumption that the gluons couple to the constituent quarks of the meson, 
and assume that the direct coupling of gluons to the V vertex is negligible. In 
general, this is not the the vertex must result from a non-perturbative 

resummation of diagrams involving quarks and gluons, to which perturbative 
gluons can in general couple directly. In fact, in all generality, such diagrams are 
needed to obtain a gauge invariant amplitude |7|, ||. However, as we shall see, 
these terms do not contribute in the high-energy limit. 

Hence we are left with the description of the qqV vertex, and the gluons 
will couple to the quarks emerging from this vertex. In general, this vertex can 
depend on the 4-momenta respectively of the vector meson, of the quark, and of 
the antiquark, V = 2v, v + 1 and v — I, as well as on 7^. We know experimentally 
that vector-meson electroproduction conserves s-channel helicity, hence the only 
possible tensor structure is 7^. We are thus left with a vertex: 

r„ = mi, (i) 

We shall assume that the vertex function $ can depend only on the relative 4- 
momentum /, and that the dependence on the meson mass can be restated as a 
dependence on a Fermi momentum scale pp. 

This vertex function can be formally related || to a light-cone wavefunction, 
but we shall not make use of such a wavefunction in the meson case. Indeed, we 
find that the analytic structure of the propagator which is usually included in the 
definition of the wavefunction plays a crucial role in the reproduction of the ratio 
vl/°~t- 111 the following, we shall make no further assumption concerning the 
quarks, except that we can treat them within perturbation theory. In particular, 
we shall not treat them a priori as on-shell particles. 
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In the following, we shall assume that the function <3> can be approximated 
for s states by a falling gaussian, multiplied by a normalisation constant N. 
As we do not assume that quarks are real particles, we do not have the usual 
wavefunction normalisation integral on $. However, the same vertex controls the 
decay V — > e + e~ , and hence we shall be able to determine N so that it correctly 
reproduces the latter rate. 



3 Calculation of the amplitude 

3.1 Kinematics 

The vectors of the problem are defined in Fig. 1. The photon has 4-momentum q 
and polarisation e, with q.q = —Q 2 . The vector meson has momentum V = q + A 
and polarisation e, with t = A 2 , V 2 = m v , hence A.q = (m v — t + Q 2 )/2. The 
quarks composing the meson are written as v + 1 and — v + / with v = V/2. 




Figure 1: The two diagrams accounting for the transition ^y*p — > Vp. The 
dashed line represents the cut which puts the intermediate state on-shell. 



As we do not want to look in detail at the diffracted proton, we follow the 
usual treatment |J: the momentum of the quark is p, and we assume that we 
can neglect its mass and put it on-shell p.p = 0. The fact that the proton quark 
remains on-shell after the scattering implies p. A = t/2. We shall be working in 
the high-u> 2 limit, and we write p.q ~ (p + q) 2 /2 = w 2 /2. We define w ~ w/3 at 
the quark level, but this will not make any difference, as the final answer, in the 
large w limit, will be independent of w. 
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3.2 The photon-meson bubble 

Assuming that we can use perturbative quarks^] in the upper bubbles, or equiv- 
alently that the non-perturbative effects can be reabsorbed into the vertex func- 
tion, the upper bubbles of the graphs are described by the following traces: 

Tf = Tr{$(/) 7 .e[ 7 .(t; + /)+m (? ]7 /3 [7.(g-t; + / + fc) + m (? ]7 a 

x h-(q-v + l) +m q ]nf.e[y.{-v + l) + m q }} (2) 

Tf = Tr{$(/) 7 .e[7.(t; -/) +m g ]7 a [ 7 .(w - k - I) + m 9 ] 7 .e 

x [^.(v-q-k-l)+m q ]j p [^.(-v-l)+m q ]} (3) 

One of the quark lines connecting to the photon is off-shell and has different 
expressions in each diagram. Its propagator in the first graph is 

2 

P x = ( q - v + l).(q - v + I) - m 2 q = (l.l + 2l.q - 21. v - 2q.v + ^-m 2 q -Q 2 ) (4) 

whereas in the second graph it reads: 

P 2 = (-v + l + k).(-v + l + k)-m 2 q 



m 



2 

V 2\ 



= (jfe.jfc + 2k.l - 2k.v + l.l - 21. v + - m 2 q ) (5) 

As we expect the amplitude to be w-independent, we shall be calculating the 
discontinuity of the amplitude, and the two intermediate quark propagators are 
put on-shell, which gives us the relations: 

k.q = A.k + A.g — k.k — 2k.l — 2l.q (6) 
l.q = -A.l + l.l + m 2 v /4:-m 2 q (7) 

One of the quarks constituting the meson is in general not cut, hence we have 
another propagator 

P 5 = (v + I) 2 -m 2 q = 2(1.1 + m 2 v /A - m\) (8) 
The sum of the two cut diagrams will then give: 

T a P _ { Tf , Tf, (2n) 2 5(P 3 )5(P±) (Q , 
~ [ Pi P 2 1 P 5 [ ' 



^One could easily extend this formalism and include non-perturbative effects present in the 
quark propagator, by replacing the quark mass vr? q by the parameters that enters in the Hilbert 
transform of the propagator, and by integrating over the Hilbert density afterwards. Such a 
complication is not needed to describe the data. 
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3.3 Current conservation 



One of the consequences of gauge invariance is that the contraction of the up- 
per bubble with the momenta of external lines is zero. In our case, we obtain 
indeed k a T a ^ = 0, and T af3 (e — ► g) = 0. However, we do not obtain exact 
current conservation for the off-shell quark emerging from the meson vertex: 
(k - A)pT a/3 ~ l.v ^ 0. In the photon case, one can recover explicit current 
conservation at the j3 vertex by crossing both sides of the diagram |l(J . This can 
be accomplished by the substitution T — > \\T(l) + T(l — > A — k — I)). In our 
case, however, because of the presence of the vertex function, such substitution 
is not possible, and the current conservation that could be checked explicitly at 
the f3 vertex in the photon case is in general absent. 

We can however impose that current conservation be satisfied, which amounts 
to putting the other quark making up the vector meson on-shell. The amplitude 
then splits into 2 physical processes: ^gg — > qq and qq — > V. This of course 
leads to current conservation and ensures gauge invariance. This is the approach 
usually followed in a wave function formalism. It has the drawback of relying 
heavily on the existence of a mass shell for the quarks, and as we shall see does 
not lead to observed results. 

However, we can pursue another route, by keeping the gauge dependence of 
the propagators, and checking explicitly that it cancels out in the large w 2 limit. 
As we shall see, the leading terms come from terms in which the momenta of the 
upper bubble get contracted with momenta of the lower quark lines. This means 
that terms proportional to g^ in the propagators are enhanced by a factor w 2 
with respect to terms containing k^ or k u , as the gluon momentum k will turn 
out to be transverse. Hence the large-w 2 limit is gauge invariant. We shall show 
explicitly that the infrared cancellation which is expected from gauge invariance 
occurs explicitly, albeit after integration over the photon-meson bubble. 



3.4 The full amplitude 

In order to obtain the full amplitude, we need to add the contribution of the 
proton lines, and the gluon propagators. We represent the proton by a constituent 
model, introducing the two following forms factors |5], 0]: 

c n _ a2\ „ (3-53 - 2.79t) 

5 ^- A )~ ( 3 .53-t)(l-t/0.71)3 (10) 

when both gluons hit the same quark line (this is a fit to the measured Dirac 
elastic form factor), and 

£ 2 (k, k - A) = £\(k 2 + (k - A) 2 - c k ■ (k - A)) (11) 



if the gluons hit different quark lines, with c ~ 1 |TT 



The prescription to go from the quark-level process to the proton-level process 
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is then to multiply the amplitude by the form factor T(k, A) = 3(£i — £2)- The 
leading contribution of the lower quark line being 4p a p lS , we obtain the following 
expression for the amplitude: 
2 

A = -(4na s ) 2 g elm e Q 

r dH r d 4 k WI AN 1 „. 

[4( W )(<T + \££W>> + 

x k'(k-Ay {U) 

where A is the gauge parameter and eQg e i m = eQ^4na e i m the electromagnetic 
coupling of the different vector mesons: tQ = ^= for the p, —1/3 for the and 
2/3 for the J/ip. 

4 The high-energy limit 

To take the large-w 2 limit, we shall use Sudakov variables, rewrite the vectors of 
the problem in terms of p and q, and define the transverse direction, noted with 
a "t" subscript, as being orthogonal to both p and q. 

We first write the expression for the polarisation vectors. The transverse 
photon polarisation is simply a unit vector in the transverse plane, whereas the 
transverse polarisation of the vector is slightly more complicated, as the trans- 
verse plane has been defined w.r.t. q. Solving e.e = —1, e.v = 0, e.e = 1 and 
e.q = 0, we obtain: 

w z Q 

e T = e t (14) 

A t -m 2 v + Q 2 -t t + w 2 
e L = + o p+ 2^ 15 ) 

e T = e t--^2~P ( 16 ) 

Taking into account the on-shell conditions, we obtain the following expan- 
sions for the measures and the 4-vectors: 

d 4 A 5 + ((p - A) 2 ) 5 + ((q + A) 2 - m 2 v ) = JL d 2 A t 
a a . m 2 v + Q 2 -t t 

A = A t + - p+^q (17) 

dH S + ((v -l).(v-l)- m\) = 8{1 \ p) d(3d 2 l t 9(1 - (3) 
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I — Q 

2 2 1 



(2A t -l t ).l t + (3(((3-l)Q 2 



my 



-t)-m v + {i 2 q 



2(P- 1) w 2 



P 



(18) 



With fig 



2m q . 



d*k 5 + ((q-v + l + k) 2 - m\) 5 + ((p - k) 2 ) = 1 ^ d 2 k t 6(1 + (3) 



k = h& < 19 > 
2[(1 - /3 2 )Q 2 - t + Afg ~ - 2A f )] - Mj + 2Z t - 2A t )(l - g) 

2(l-/5 2 )w 2 P 



The propagators then become: 



P 5 = 2Z 2 + ^-^ 

P 

1.2 _ ?t_ 
4 

(A; - A) 2 = -A t .k t + ^ + t (20) 
In the following, we shall write the integrated amplitude as: 

= ^(4tt«5) 2 gei m e Q 
d 2 k r dH 

X 



r ark r 
J J2~kVJ 



(2tt) 2 J (2tt) 4 
x F(k,A) -±=<S>(1) A (T , L) (21) 

The leading terms in w 2 for the two amplitudes -A(t,l) then become A T = 
N T /D and A L = N L /D with 

£> = (2(3A t .l t -l t .l t + ^ 2 q + ((3 2 -l)m 2 v -t(3 2 ) 

x ((1 - f3 2 )Q 2 + fJ , 2 g-t-(l t + k t - 2 A t ). (k + k t )) 

x ((i-p 2 )Q 2 + (2A t -l t ).l t + f i 2 q -t) 

x (4A t ./q + /q./q-4t)/q./q (22) 

N L = 4(1 g K(/3 2 - 1) + t/3 2 - ^ + l t .(l t - 2(5 A t )] 

x k t .(2l t + kt - 2 A t ) (23) 



9 



N T = 8e tfl e tu 

x {[-A? A?(/J + 1)/? + - ZfZr(/? + 1) + *W(/3 2 + (3 + 1)] 

(fct + 2Z t -2A t ).A; t (l-/3) 
+[-/3AfA;r - Kr t (3 2 + A*/? 3 + 
x[t- t i 2 q -(l-p 2 )0 2 + h.l t -2A t .k] 
-g^[-k t .(k t + l t ){lflt - $ - (1 - P 2 )k t .l t Q 2 
+A t . ((^ - Z t .Z t )(/3 - 2)A* + + l)fc t .A« + 2/3^J t - 2A t .fc t )Z t ) 
+tk t .(A t f3-ktf3-2ltP + lt) 

+(l-(3 2 )(3A t .k t Q 2 }} (24) 




In order to proceed further, we choose two orthogonal polarisation vectors 
in the transverse case, and there are a priory four transverse amplitudes, as the 
polarisation vector can each be taken parallel or orthogonal to a fixed direction. 
However, we find that the non diagonal terms, in which the photon and the vector 
meson have orthogonal polarisations, cancel in the angular integration. Hence in 
general we have only two diagonal transverse amplitudes to consider. 

5 Properties of the amplitudes and of the cross 
section 

To pursue analytically, we shall concentrate on the case t — and A = (my + 
Q 2 )/w 2 p. We shall give numerical results at all t further on. We can then perform 
the angular integrals analytically. In this case, the remaining two transverse 
amplitudes become equal after angular integration. Hence we are left with a 
transverse and longitudinal amplitudes, dAi and cIAt, which we still need to 
integrate over P, k t 2 and l t 2 . We use bold-faced letters to denote euclidian 
vectors in the transverse plane, hence k 2 = — k t 2 , etc. 

We cannot perform the k t 2 integral, as it depends on the proton form-factor. 
For the l t 2 integral, we change variable to Z 2 as this is the variable entering the 
vertex function. 



As we are going to discuss the real part and the imaginary part of the expression, 
we express the pole explicitly: 




l t 2 = (l-/? 2 )m 2 / - / iJ + 2(/^-l)P 



(25) 



V = 2Z 2 + 




+ ie 



(26) 



Defining : 
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L 2 = (1 - p 2 )(Q 2 + m 2 v ) - 2(1 - P)V 
M 2 = 2(1 - {3 2 )m 2 v - 2(1 - /3)V 
ir 4 = (L 2 + k t 2 ) 2 -41 2 k t 2 

we obtain after angular integration: 

- 

dA T = — — — -i— , d/ 2 d/3 rfk t 2 

(2tt) 3 K 2 L 2 k t 4 2P 

{(K 2 - L 2 ) [(1 + f){L 2 - 21 t 2 ) - 4/i 2 .] - k t 2 L 2 (l + f) } 

(28) 

The positive energy condition (v — l).v > mv ^ q implies that 

I 2 < {mv ~ ^ (29) 

These expressions have the following remarkable properties: 

i) The small- k t behaviour of both amplitudes is precisely l/k t 2 . This is essential 
as it means that there is no IR singularity: the remaining divergence is handled 
by the proton form factor. We get: 

ii) One recovers the ratio Al/At of our previous work, in the limit of zero Fermi 
momentum: 

2fi 2 dl 2 d(3dk t 2 , 

dMfi = °' 0) = WVT^W <32) 

<M ^ = 0,I, = 0) = (33) 



I ci4 z 



For yU 2 = my and Q » m v , we recover 
iii) There is no divergence at the edge of the (3 integrals. The longitudinal cross 
section is zero, and the transverse cross section goes to a constant, independent 
ofQ 2 . 

dA T (0 = -l) = %^t^V dl2 dpdkt2 (34) 



1Aln , dl 2 dp dk t 2 
dA T ((3 = +1 = - — xo * x 
KH ! 2tt 3 k t 2 V 



1 (m 2 v -2Q 2 -Al 2 + fi 2 ) 



+ 



K 2 ' k t 2 (m 2 / + 2Q 2 -4/ 2 + /i 2 )J (35) 
dA L ((3 = ±1) = (36) 
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We can now perform the j3 integration. Taking into account condition 
we write I 2 = (my — /i g ) 2 /4 — A 2 . The positivity of \ 2 from (f2~5|) constrains j3 to 
be contained between the bounds: 



±2-,/A 2 + ii a my \ — 2A — a a my + m\ 
P± = — * 5 " ( 37 ) 

The lower bound of integration is always comfortably below 1 for nonzero \x q and 
finite A. The lower bound becomes smaller than -1, and is then not realised, for 
A > \jm\ — m q my. 

However, the amplitude has a pole in I 2 , or equivalently in A 2 , at 

V = o =► I 2 = ti-™* (38) 
4 

hence we obtain a contribution from the principal part integration, and a con- 
tribution from the discontinuity. The latter corresponds to a direct extension of 
our previous model: both quarks constituting the meson are on-shell, and, after 
integrating over f3, we obtain the following for the contribution to the amplitude 
at high Q 2 : 

-8m v p d dl 2 dk t 2 

dA L[dlsc) « -^qT^T (39) 

-4m&fa + 2(/i 2 + 2m 2 ,) log(f±|) d p d k t 2 
dA T{disc) - ^ ~ W ( 4 °) 

2 2 

The bounds on j3 corresponding to (p7|) are /3J = my m i iiq , far away from (3 = ±1. 
We see that the ratio dAi/dAx is still linear, but depends on the quark mass 
chosen. In the limit (3d — > 0, we of course recover our previous results for the 
ratio, and that the ratio Q/my gets multiplied by a constant, which depends on 
the quark mass. Unfortunately, this constant is always between 0.5 and 1, as 
shown in Fig. 2 for reasonable values of the (constituent) quark mass. Hence the 
discontinuity alone cannot solve the problem of the ratio <Jl/&t- 

Note that it is possible to extend somewhat the previous expressions. Al- 
though the expressions include Fermi motion, they really single out a value of I 2 
which corresponds to the on-shell condition, and the remaining f3 integral can be 
reduced to an angular integral in the meson center-of mass frame. One can make 
another model, where the quark mass is a priori not fixed, but a function of I 2 , 
chosen so that the on-shell condition ( |3~8D is satisfied. The remaining integration 
involves then the vertex function $ (I), and it cannot be done analytically. We 
have tried this route numerically, and find that the ratio remains linear, and 
can be somewhat reduced, by about a factor 2. This is not enough to repro- 
duce the measured values, and is reminiscent of the situation encountered using 
a light-cone wavefunction formalism. 
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Figure 2: The reduction factor in the ratio cl/ctt as a function of the 
quark mass 



At this point, we still have not included the contribution of the principal 
part integration. This corresponds to letting one of the quarks off-shell, and 
there is no reason to neglect it. In fact, in the case of DIS, the kinematics is 
such that this is the only contribution. Its structure is somewhat different from 
that of the discontinuity. Indeed, the value of I 2 is not fixed anymore, hence 
the integration over f3 has only the bounds ([37]). In particular, once I 2 becomes 
large and negative, the value (3 = — 1 can be realised. As we have mentioned in 
Eqs. (|35D, the transverse part retains a finite value there, whereas the longitudinal 
part goes to zero. Hence there is a narrow peak which gives a contribution only 
in the transverse case. If we include only the large A contribution (the other parts 
turn out to be negligible), one obtains: 



dA 



dl 2 dk f 2 



L(PP) 



(2vr) 3 V m y Q 2 



+ log 



(4P - k t 2 ) V 

(i + /3+) 2 Q 4 



2m 2 / (l + /? 4 
V 

V 

-4-^loe 



4P - k t - 
2V 



dA T (pp) 



-2 log( 



4P-k t 2 
4V 



dl 2 dW 



k t 4 Q 2 



(27T) 



(41) 
(42) 



We can see that the real parts behave like dA^ oc and dA? oc ^2 at high 
Q 2 , plus logarithmic corrections, whereas the imaginary parts still behave like 
dAi oc ^3, and dA T oc Because this effect is present only at large A, it 
is suppressed by the fall-off of the vertex function, and sets in only at relatively 
large Q 2 . The leading behaviour of the principal parts integrals comes from quark 
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off-shellnesses bigger than the constituent quark mass. Whereas in the massless 
case || an increase in the cross section can only come from extremely small off- 
shellnesses, we find that the dominant region in the massive case is shifted by 
the quark mass: the principal part integral cancels as long as one is very close to 
the pole, and the kinematics is such that this cancellation is not present anymore 
once the off-shellness is of the order of the quark mass. 

Hence the source of the plateau observed at HERA is the interplay between 
the real part and the imaginary part of the cross section. Our model in fact 
predicts that asymptotically the transverse and the longitudinal cross sections 
first become equal, and that ultimately the process is dominated by the transverse 
cross section. This prediction is however driven by the details of the quark 
propagators, which could be modified by confinement effects. 



6 Results 

We shall now give our numerical results, which come from a numerical study of 
the full integrands at all t values. Before doing so, though, we shall need to define 
the vertex function, and to normalise it. 

6.1 Vertex function and normalisation 

The function $ is unknown, and only its general analytical properties are well 
established ||. In general, this function can depend on the scalar products V.V, 
V.l and l.l. However, as in our case one of the quarks is on-shell, and as the 
meson is on-shell, only one of these scalar products is free. As in the case of the 
imaginary part, our results reduce to those of a wavefunction formalism, we shall 
assume that the vertex function is similar to the wavefunction of an s state. In 
the case of the p, the J /if) and the 0, we take the form: 



$ (/) = AT e 2p f (43) 

where L 2 is the quark 3-momentum in the meson rest frame, equal to L 2 = 
( — ) 2 — l.l, where the Fermi momentum pp is 0.3 GeV in the p and <p cases, 



and 0.6 GeV in the J/ip case ||12|| , and where the quark masses are taken to be 
r n u ,d = 0.3 GeV, m s = 0.45 GeV and m c = 1.5 GeV. 

We can now fix the normalisation constant iV by requiring that the decay rate 
T — > e + e~ be reproduced through the process depicted in Fig. 3. 

The total decay rate, T, for V — > f\ji is given by: 



I -j^\ M \ 25{A){ y - * - w* 2 - m f 2 ^ 2 - ^ 



2m v J (2 
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-v/2+1 




v/2 + 1 



v/2+1 




v/2+1 



Figure 3: Diagram renormalising the vertex 



with M the invariant amplitude corresponding to the diagram of Fig. 3: 



d 4 l 

(27 

x (T 2 ) 

T 3 ^ m 2 e g (^)$(0 



x 



(2tt)< 



1 



(45) 



With 



Tr 
T 3 



Tr[y.e(>y.(-v + I) + m q )Y(l-(v + I) + m q )\ 
[(-v + I) 2 - m 2 ] [(v + I) 2 - m 2 } 

Tr[Y(l-f2 - m/)7 M (7./i - m f )\ 

Tr[7^(7.(-w + l) + m q )^ .e(7 .(v + I) + m q )\ 
[(-v + I) 2 - m 2 ][(v + I) 2 - m 2 ] 



(46) 
(47) 
(48) 



Putting the final leptons on-shell, neglecting their mass, working in the center- 
of-mass frame of the meson, where I 2 = Iq — L 2 , and performing the angular 
integrals, leads to: 



r = 



X 



1 ttV 2 9 



2m v (2tt) 2 
L 2 dl dL 



4 2 

9elm 



J_ 128, 

\y 2 ) V g , 



24 3 

[3(4m 2 + m 2 v ) + 4L 2 - 121 2 ] $ (/) 
(2tt) 3 (4Z 2 , + 4l m v - 4L 2 - 4m 2 + m\ + ie) 
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The integration over l can be carried out by residues. We then get: 

1 2 



2tt 



9 



2my 3 

L 2 dL 



3m v 



3(2tt) 



(2mt) 



^2L 2 + 3m 2 g ) $ (Z) 



L 2 + m2 (L 2 + m2 



(50) 



The remaining integral over L still has a pole at the quark propagator. We 
must again keep both the discontinuity and the principal part of the integral. We 
in fact see that the phase of the amplitude becomes pure imaginary again once 
the vertex gets properly normalised. Introducing the function we have to deal 



with the pole in the denominator for L = y — m 2 . 
We introduce the notations: 



z 

Zd = 

/(*) 



2ILI 



mi 



(4 



z 2 + K 



(2z* + 3^ q ) $(z 2 /4) 
(z - z d ) 



(51) 



with $(z /4) = N e F t then the integration over L gives: 



4 «ezm eg 2 
9m v n 



P 



f(z)dz 



1 2 



Zd 



+ (iTT) f(z d ) 



(52) 



It is worth pointing out that if we neglect the momentum Z in the quarks loop, 



e.g. for m q 



2 ' 



we recover the formula |J for the p meson: 



r 



jm ^_ 2 2 
o a elm e Q 
my 3 v 



(53) 



We give in the Table the values of the decay rates [I|j] which we have fitted to, 
and the corresponding values of the normalisation N. 



Meson 


r(^ - 


-> e + e~) (keV) 


|Af 


p F (GeV) 


m q (GeV) 


P 




6.77 


61.71 


0.3 


0.3 






1.37 


72.01 


0.3 


0.45 






5.26 


44.91 


0.6 


1.5 



Meson decay rates, vertex normalisation constants, Fermi momenta and 
quark masses. 
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★ HI 

□ ZEUS 
o KMC 
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Q 2 (GeV 2 



Q 2 (GeV 2 ) 



Figure 4: (a) Cross sections for p as functions of Q 2 , compared with data 
from ZEUS § and HI § at < w >«100 GeV, (b) Ratio of cross sections 
as functions of Q 2 at < w >~100 GeV, compared with data from HI and 
Zeus §, II, and from NMC p]. 



6.2 Q 2 and mass dependence of the cross section 

The Q 2 behaviour of the total cross section, as well as the mass dependence, 
hardly get affected by the addition of Fermi momentum. This is because at low 
Q 2 , the transverse amplitude is not modified by the addition of the real part, 
whereas at high Q 2 the cross section is still dominated by the longitudinal part 
at present values of Q 2 . We show in Fig. 4. a the result of our model together 
with data from HERA. We have not corrected either for an eventual mass or Q 2 
dependence of the Regge factor. Also, as our model does not include a prediction 
of the energy dependence, we have compared with 1994 data. New data are more 
precise, and allow one to observe that energy dependence directly from HERA 
data. The figure (plotted for data at w 2 ~ 100 GeV 2 ) shows however that the 
general features of the data are well reproduced: we see that a constant Regge 
factor (equal to 5.1) reproduces the data well for all mesons at HERA. This can 
be checked by considering figure 4.b which shows the ratio of cross sections, which 
goes well through the absolute predictions of our model. As in our previous model 
we see that there seems to be a deficit of 4> mesons w.r.t. p mesons at NMC, 
which we are unable to account for. Comparison with lower-energy data may lead 
to an estimate of as as well as of the intercept. These estimates are identical to 
those of although the value of as should be multiplied by a factor v^3 which 
we overlooked in our previous work @. It is again possible to go through the 
photoproduction point, although admittedly our model should not work for such 
low values of Q 2 . 

We have explained in great detail at t = that ratio of the longitudinal and 
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Figure 5: (a) Scaled behaviour of the real and imaginary parts of the 
transverse and longitudinal amplitudes, (b) Ratio of the longitudinal and 
transverse parts of the p cross section as functions of Q 2 



transverse amplitudes has a plateau at large Q 2 . It remains to be seen whether 
this feature is maintained for the t-integrated cross section. We show in Figure 5. a 
that this is indeed the case. Furthermore, as demonstrated in Fig. 5.b, for the 
values of Fermi momentum and for the quark masses shown in the Table, one 
falls right on top of the data. There is of course some uncertainty linked to these 
values, but a change of 50% in either leads to a change of about 20% in the ratio 
a L /a T . 

6.3 Photoproduction, Regge factor and t dependence of 
the cross section 

As we briefly mentioned, our model works in the deep non-perturbative region 
of photoproduction. It may be worth mentioning at this point that even at 
large values of Q 2 the gluon off-shellnesses are not large, and as we have seen 
quark off-shellnesses must be near m q in the real part. Hence the situation is not 
dramatically different in the case of photoproduction, but nevertheless, it still 
comes as a surprise that the model applies in this region. This may hint at a 
lower-Q 2 generalisation of the factorisation theorem. 

We show in Fig. 6. a the p photoproduction cross section da/dt and compare it 
with our model. We see that we obtain excellent agreement, although our curve is 
not an exponential. In Fig. 6.b, we show the Q 2 dependence of the b slopes from 
our model: the experimental points correspond to a fit to Ne bt , and our curve 
corresponds to l/< t >, which would be the same were the curve an exponential. 
We see that despite the curvature of our curves, the agreement is quite good. The 
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Figure 6: (a) Differential cross section da/dt for elastic p° photoproduc- 
tion, compared with HERA data Q], (b) i-slope of the photoproduction 
differential elastic cross sections, compared with HERA data 



t-dependence of other mesonic cross sections is predicted in Fig. 7. Comparison 
with preliminary "public" data from ZEUS ]15| would indicate good agreement. 



We trust that the forthcoming ZEUS published analysis will confirm this result. 
Hence, if we believe this model is applicable to photoproduction, we see there is 
little room either for a BFKL large t enhancement or for a pomeron slope. 

6.4 Predictions for p' and ty' 

Given the success of our model in reproducing the lowest mass vector mesons, 
we can easily extend it to study the 2s excited states. The first ingredient is the 
vertex function, which we take as: 

3 



* 2 , = AV57i(- 5 + ^jexp (54) 



with L as defined below Eq. (|43|), and the values of pp and m q the same as in 
the Is case. The normalisation constant N should again be determined from the 
leptonic decay of the vector meson. The latter is available only in the ip' case 
JTJJ], hence we can make predictions only in this case: 7(/0' — > e + e _ ) = 5.26 keV 
leads to N = 35.36. The ratio of the ip' elastic production to that of the J/ip, 
at low Q 2 < 0.01 GeV 2 is calculated to be 0.165, in excellent agreement with the 



HI measurement fl6|| , as can be seen from Fig. 8. a. 

In the p' case, we can calculate from the elastic production and make 
predictions for the leptonic width: from the preliminary "public" results of the 
HI collaboration [jlTH , which we were unfortunately asked not to quote |Q, we 
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1*1 (GeV 2 ) |i| (GeV 2 ) 

Figure 7: Ratio of the differential photoproduction cross sections for (a) <f> 
and p mesons and (b) for J ftp and p mesons. 




Figure 8: (a) Prediction for the ratio of the ip' production cross section to 
that of the J/ip, compared with HI data (b) prediction for the p' cross 
section based on the consideration of HI preliminary "public" data 17], 
leading to the prediction of the leptonic width (|55|). 
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obtain 

T(p' -> e + e~) = 1.1 ± 0.34 keV (55) 

The resulting predictions are shown in Fig. 8.b. We see that although the nor- 
malisation of this figure is not predicted, the ratio is predicted to increase sharply 
with Q 2 and hence an eventual final measurement at HERA should find a sub- 
stantially higher cross section than the one at fixed target experiments. 

7 Conclusion 

We have presented here a model which reproduces all the features (except the w 2 
dependence) of elastic vector meson production as observed at HERA. Although 
the model is built to work at large values of Q 2 , it extends to the nonperturbative 
photoproduction region. The main result is that it is possible to reproduce the 
ratio OlIot- The plateau observed experimentally comes from the interplay 
between the on-shell and the off-shell quark contributions, which have different 
asymptotic behaviours. Hence it is essential to allow quarks to be virtual, as in 
the standard DIS case. 

It is worth pointing out that despite its successes, our model has not been 
tuned to reproduce the data: the values of the Fermi momentum, the form of 
the wavefunction and the quark masses are given in the Table, and are only 
reasonable guesses. Similarly, the proton form factor could be modified, as it is 
known only in the IR limit. We find remarkable that educated guesses lead to 
such a good agreement with data. 

Hence it is essential to allow quarks to be off-shell not only to reproduce 
the behaviour of the transverse part, but also because the relation to structure 
functions is only possible then. Indeed, the logQ 2 terms of Eq. ( fill) are the 
same as those of the upper quark loop in DIS. The use of wavefunctions can miss 
such contributions, as they come from a part of the amplitude which is usually 
embedded in the wavefunction. 



Note 

A previous version of this paper showed preliminary HI and ZEUS data scanned 
from "public" documents available from the web and comparing them with our 
model. Although we found beautiful agreement, the HI collaboration demanded 
that we remove these from the published version, and the ZEUS collaboration 
indicated that they as well do not want preliminary results included for compari- 
son in theoretical papers. We have thus removed these data, and deeply regret to 
have taken the scientific liberty of explaining experimental results publicly avail- 
able before they are published. We want to indicate that that policy should be 
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spelled out explicitly in the papers and in the webpages available, as comparing 
results with (even unpublished) data is -and should be- a common trend among 
theorists. 
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